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Abstract 

New method for dissociation of diatomic molecule based on nonperiodic exci- 
tation generated by energy-feedback control mechanism is proposed. The energy- 
feedback control uses frequency-energy (FE) relation of the natural oscillations to 
fulfill the resonance conditions at any time of excitation. Efficiency of the proposed 
method is demonstrated by the problem of dissociation of hydrogen fluoride (HF) 
molecule. It is shown that new method is more efficient then methods based on 
constant frequency and linear chirping excitation. 
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1 Introduction 

During last decade a growing interest has been observed in the control problems for 

molecular systems in classical and quantum formulation 0, @, U- One of the simplest 

problems of that type is dissociation problem for diatomic molecule || f|, ||, [| |7|, ||. It 

often serves as a benchmark for comparison of classical and quantum calculations. In the 

paper ||] possibilities of dissociation a molecule by monochromatic (single frequency) laser 

field have been explored for the case of hydrogen fluoride (HF) molecule using Chirikov's 

resonance overlap criterion. In || U the case of two-frequency (two-laser) control field was 

investigated. It was shown that intensity of two-frequency field required for dissociation 

can be reduced compared to one-frequency case. In |J the further reduction of the control 
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field intensity has been demonstrated by means of chirping (frequency modulation) the 
laser frequency with constant chirping rate. 

New possibilities for changing of physical and chemical properties are provided by 
using feedback. In || methods based on geometric control theory (inverse control) were 
proposed for molecular systems. In the papers P, [II]] a general method for investigation 
of physical systems by means of feedback controlling actions was proposed. It was shown 
in || by example of the problem of escape from a potential well that feedback allows to 
reduce control intensity required for overcoming a potential barrier by several orders of 
magnitude. It appears that using speed-gradient method fly], |l2j for design of feedback 



control action allows to achieve desired level of energy by excitation of intensity inversely 
proportional to the system dissipation degree. Particularly, for conservative systems it 
allows to reach any energy level by (ideally) arbitrarily small control action. 

In this paper the algorithm of the feedback frequency chirping of the controlling field 
frequency is proposed based on dependence of the natural frequency of the nonlinear 
molecule oscillations on its energy, further referred to as frequency- energy characteristics 
or frequency- energy curve (FE-curve). On-line implementation of the algorithm requires 
measurement of molecule energy only. This can be easily generalized for an ensemble of 
molecules, replacing exact energy of the separate molecule in the control criterion by the 
average energy per single molecule. However, for control algorithm design the knowledge 
of the frequency-energy relation is required. In this paper a classical Morse potential is 
considered, which allows obtaining this relation in a simple analytical form. For other 
potentials, where the analytical solution is too complicated or impossible, the frequency- 
energy relation can be obtained numerically. Efficiency of the proposed algorithm for the 
case of HF molecule is illustrated by computer simulation. It is shown that the time of 
dissociation can be reduced significantly when using proper feedback excitation. 

2 Basic equations 

Let us consider a diatomic molecular system under the action of the external laser field. 
Dynamics of such a system can be described by the following controlled Hamiltonian || [3|] 

2 

H = ^ + Il(r)-n(r)u(t), (1) 

where m is reduced mass, p is momentum, IT(r) is potential of interatomic interaction, 
/i(r) is dipole moment of the molecule, u(t) is intensity of external field. The value u(t) 
serves as control variable. E.g. for monochromatic control field u(t) = Ecos(ut), where 
E and uj are strength and frequency of the external laser field. Substitution of ([!]) into 
Hamilton equations 

dr dH dp dH 
dt dp ' dt dr 
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yields the following equation of molecular motion 

mf = -Il'(r) + n'(r)u(t) . (3) 
For the case of harmonic excitation equations of motion is as follows: 

mr = — n'(r) + Efi'(r) cos(ut) . (4) 
Let us use Morse form of the interatomic potential 

n(r) = D (l - e - a{r - a) ) 2 — D = D ( e - 2a ^ - 2e- Q(r - a) ) , (5) 



where D is the bond energy, a is the equilibrium interatomic distance. The corresponding 
force is 

f(r) = -n'(r) = 2aD ( e - 2a(r " a) - e - a(r " a) ) . (6) 
The dipole moment can be represented in the form || |13[ 

/i(r) = Are"^ 4 //(r) = A (l - 4£r 4 ) e"^ 4 , (7) 

where A and ^ are the known constants. Thus the equation of motion |] reads 

mr = 2aD [e- 2a{r - a) - e^^) + EA (l — 4£r 4 ) e^ r4 u{t) . (8) 

For sufficiently small £ the simplified expression for dipol potential can be used: 

H(r) = Ar, fj/(r) = A, (9) 

In the vicinity of equilibrium r pa a equation (|]) with harmonic excitation can be reduced 
to the equation of linear forced oscillations 

mr = Cr + Efi'(a)cos(ujt) , C = II" (a) = 2a 2 D , (10) 

where C is linear stiffness of the bond. The natural frequency of the linear system is 



Q d 4 f yfc/m = a^2D/m . (11) 

3 Control of the dissociation process 

In a nonlinear system the resonance frequency is function of the amplitude of oscillations, 
or, in other words, it depends on the energy of the oscillator. Thus the faster dissociation 
can be achieved if the excitation frequency is changing (decreasing in our case) while 
the energy increases. The simplest control algorithm of this kind can be done by linear 
chirping: 

uj = n -et (12) 



where e is a constant (chirping rate), which characterizes the speed of the frequency 
decrease. In this case equation of motion ([J) takes the form 

mr = -W (r) + Efi'(r) cos <p , tp = u(t) = O - et . (13) 

(Note that in this case the representation tp — ut is not valid since the frequency u is not 
constant. 

More effective dissociation can be achieved by taking advantage of the feedback control. 
In this case the frequency of excitation can be changed according to the oscillator energy 
so that at any amplitude of oscillations the excitation will act at the resonance frequency. 
Let f)(W) be the natural frequency of the diatomic molecular oscillator at the specified 
energy W. We call the function f)(W) frequency- energy curve of FE-curve. It contains 
important information about dynamics of the molecular motion. The key idea of our 
approach is that if the curve fl(W) is known and the energy W = W(t) is observable then 
control algorithm 

mr = -n'(r) + E/j'(r) cos tp , (p = (14) 

will give the desired excitation at the resonance frequency at any time instant. 

The relation ^(W 7 ) in principle can be calculated for any known potential U(r) using 
the integral of energy, which is realized when the external force is absent 

mf 2 + Tl(r) = W => T = V2m I — (15) 



2 ;/ yfW-U(r 

where W is a constant value of energy, T is period of the oscillations, r\ and r<i are the 
minimum and the maximum possible values of radius r for the specified potential energy 
(the solutions of the equation U(r) = W). If period T is known then frequency can 
be obtained as Q = 2n/T. Obviously integral (O) can not be calculated analytically 
for arbitrary potential and the desired relation Q(W) should be evaluated by numeric 
integration. 

Luckily for the Morse potential (^) the integral can be evaluated analytically in the 
following way. Substitution £ = £(r) = e a(r_a - ) transforms the integral fllBT) in the case of 
Morse potential to the simple form which can be easily integrated 

C2 / 

d( = /2m 

K C 2 + 2C _ i a V W 

The above formula is valid if W < 0. Otherwise, the motion of the system is not peri- 
odic, since the molecule dissociates. Equation (|T6| ) provides the following simple relation 
between energy and natural frequency of Morse oscillator 




15 ~ vr) = \ fi~ 



(17) 
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where Qq is defined by (IT]), T = 2tc/Q . The control law ([14]) therefore takes the form 

^ = n(w) = n 0] j . (is) 

where W(t) is the current value of the molecule energy. 

4 Calculations 

Efficiency of the proposed control algorithm was confirmed by computer simulations for 
dissociation of the hydrogen fluoride (HF) molecule that has become a benchmark problem 
for controlled dissociation || || The parameters of the model [8], corresponding to HF 
molecule i, i are as follows: m = 1732, D = 0.2101, a = 1.22, a = 1.75, A = 0.4541, 
£ = 0.0064, E = 0.1. All quantities are given in atomic units. The initial condition for 
excitation is molecule equilibrium at the bottom of the potential well: r = a, r = 0. 
The equation was integrated numerically by central differences method with time step 
dt = Tq/ 160 (integration with smaller time steps yields practically the same results. 




Figure 1: Excitation at constant frequency: a) resonance frequency uj = Q , b) optimal 
frequency u = 0.8A1Q . Energy is measured in the units of the bond energy D, time is in 
the units of T (natural period of linear oscillations). 



The result of excitation at the resonance frequency is shown in Fig. la. The graph 
shows that the motion is quasiperiodic and maximum value of energy, which can be 
achieved by this method, is — 0.68-D. Since the natural frequency of the molecule is 
decreasing when the energy is growing, more efficient excitation can be achieved at lower 
frequencies. It is even possible to dissociate the molecule if the excitation frequency is close 
to 0.841f2o — see Fig. 2a. The figure shows that at t ~ 13.5To the energy crosses the zero 
level, and at t = 14.5T the energy becomes a positive constant. This practically means 
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that the potential energy is negligible and atoms are moving in the opposite directions 
with constant velocities (total dissociation). However this result is not structurally stable: 
small variation in the excitation frequency can noticeably increase the dissociation time 
and even break the dissociation. 




Figure 2: Excitation with linear frequency chirping: a) e — 0.0157 ^r, b) e = 0.059 
(the chirping rate, which provides the fastest dissociation). Energy and time are measured 
in units of D and To. 

Fig. 2a,b correspond to linear frequency chirping (|12"D at two different chirping rates. 
The dissociation is also achieved, but it is slower then in the previous case and the results 
are very sensitive to the variations of e. Additional flexibility to the chirping process can 
added by using the relation to = lu q — et, where uq is a constant less then Q . Optimal 
selection of both parameters is a pretty complicated task. In more detail this problem is 
investigated in [||. 

In this paper we will concentrate on the feedback chirping, based on the frequency- 
energy curve (|TT| ) . Fig. 3a shows the energy as a function of time for the feedback con- 
trol ([T8|) . At the initial stage the energy is growing faster then in other cases, but after 
the level about W = —0.5D is reached, the efficiency of the excitation decreases so that 
dissociation is not realized. Probably this phenomenon is caused by wrong phasing which 
prevents the molecule dissociation although the resonance condition is all the time fulfilled. 
Indeed, the sign of the multiplier /i'( r ) i n equation (||) due to formula (0) is changing for 
high energies, so that the excitation starts acting in the opposite phase. To overcome this 
difficulty the following algorithm can be used: when the energy of the molecule reaches 
some certain value W* then the sign of the excitation force is permanently changing to 
the opposite. The result of such an excitation is shown in Fig. 3b. The dissociation is 
achieved at t = 10T , which is faster then in all previous cases. 
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Figure 3: Excitation with feedback control: a) using the frequency-energy relation; 
b) same, but with phase switch at W = —0.2D. Energy and time are measured in D 
and T - 

Thus for the selected level of control (E = 0.1 a.u.) the best result is obtained 
for feedback excitation. However excitation with shifted constant frequency gives the 
competitive result. The situation changes if we consider lower level of excitation. If 
this level is sufficiently small then excitation at any constant frequency can not provide 
the dissociation. On the contrary, excitation with changing frequency can in principle 
give dissociation for any arbitrarily small level of control. For illustration let us consider 
the 5 times smaller control then in previous calculations: E = 0.02 a.u. In this case 
dissociation at a constant frequency is not possible. Results of linear chirping and feedback 
control excitation are given in Fig. 4. For simplicity these calculations use approximation 
//(r) fi'(a). The plots show that both methods provide dissociation. However, linear 
chirping achieves dissociation at t — 89T , while energy-feedback control achieves it at 
t = 34T . Thus the feedback control provides more then two times faster dissociation 
then the best variant of the linear chirping. 

Although the proposed algorithm requires measurement of molecule energy only, its 
on-line implementation requires changing controlling action at the time intervals less then 
10~ 12 s. This can be achieved with modern ultrafast lasers. The approach can be easily 
generalized to control of an ensemble of molecules, replacing exact energy of the separate 
molecule in the control objective by the average energy per single molecule. 
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